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ABSTRACT 

We employed published rotation periods of Kepler field stars to test whether stars 
hosting planets tend to rotate more slowly than stars without known planets. Spec¬ 
troscopic vsini observations of nearby stars with planets have indicated that they tend 
to have smaller visni values. We employ data for Kepler Objects of Interest (KOIs) 
from the first 16 quarters of its original mission; stellar parameters are based on the 
analysis of the first 17 quarters. We confirm that KOI stars rotate more slowly with 
much greater confidence than we had previously found for nearby stars with planets. 
Furthermore, we find that stars with planets of all types rotate more slowly, not just 
stars with giant planets. 

1 INTRODUCTION 2 PREPARATION OF SAMPLES 


Stars with planets (SWPs) have been shown to be more 
metal-rich as a group when compared to similar stars with¬ 
out detected planets (non-SWPs) for nearby stars, found 
with the Doppler method (Valenti and Fischer 20051, and 
for stars in the Kepler field (Ji and Fischer 20151, found 


with the photometric transit method. SWPs and non-SWPs 


|Alves et al.|2010 Delgado Mena et al.|2015|. 


have also been reported to differ in vsini (Gonzalez 2008 
|Gonzalez et al. 2010 Takeda et al. 2010 Gonzalez 20111, 


but others have failed to confirm this (|Israelian et al.||2009| 


Bouvier (20081 suggested that slow rotation among 


SWPs caused by early star-disk interactions caused accel¬ 
erated Li destruction in them. Recently, additional evidence 
for lower Li abundances among SWPs has been presented 
by Figueira et al. (20141 and (Gonzalez (20151, which, if the 
evidence for slower rotation among SWPs is also confirmed, 
strengthens the Li abundance, stellar rotation and planet 
links. Before Bouvier’s theory would be accepted as the best 
explanation, however, it would have to be shown that the 
spin-down of stellar hosts occurr very early. 


The purpose of the present study is to test whether 
SWPs rotate more slowly non-SWPs using Kepler data and 


a method of analysis similar to that described in (Gonza- 
lez|2015 I. The paper is organized as follows. In Section 2 we 


describe the preparation of the SWP and non-SWP compar¬ 
ison samples. In Section 3 we compare the measured rotation 
periods and simulated vsini values of SWPs and non-SWPs. 
We discuss the result sin Section 4 and present our conclu¬ 
sions in Section 5. 


In order to compare the rotation periods of non-SWPs and 
SWPs in the Kepler field, we employed the results of the 
McQuillan et al. | (|2014[ ) and |Mazeh et al.| ( |2015| ) studies. [Mc-| 
Quillan et al. (2014I determined rotation periods for 34,030 


main sequence targets in the Kepler field using observations 
from the Q3-Q14. Their dataset also excludes known eclips¬ 
ing binaries and KOIs that had been identified by Septem¬ 


ber 2013. In a follow-up study, Mazeh et al. (2015) applied 


the same analysis method as McQuillan et al. to the Kepler 
KOIs and published 3,356 rotation period values for the sub¬ 
set with detected periods. Hereafter, we will refer to these 
two datasets as McQ14 and Mazl5, respectively. We also 
downloaded (on February 3, 2015) the latest list of KOIs 
from Q1-Q16 (hereafter, KOI16) and the stellar parameter 
determinations from Q1-Q17 ’Kepler Stellar Table’ (here¬ 
after, S17) from the Kepler archives]^] KOI16 includes 7,349 
stars, and S17 includes 199,244 stars. 

Since one of our goals is to compare stellar rotation 
periods in the Kepler field stars to visini values measured 
in nearby stars by Gonzalez (20111, we must restrict the 


stellar parameters of the Kepler samples to the same range 
as the nearby stars sample. Thus, only stars with effective 
temperature (T e ff) between 5500 and 6500 K were retained 
in our four Kepler datasets; exclusion of stars cooler than 
5500 K also reduces the chances of including giants in the 
sample. In addition, only stars with log g values greater 
than 3.5 were retained with the express prupose of excluding 
giants. We removed stars with assumed solar parameters and 
stars lacking T e ff or log g values in the S17 dataset. This 
resulted in 107,633 stars. The stellar parameters from S17 
replace the stellar parameters listed in the other datasets. 


http: //exoplanetarchive.ipac.caltech.edu 
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2 G. Gonzalez 


Typical quoted uncertainties in T e ff, log g, and [Fe/H] are 
^170 K, io. 2 7 dex, and IqJs dex, respectively. 

After removing any known KOIs from the McQ14 
dataset (using the ’nkoi’ column in the S17 dataset), it con¬ 
tained 15,467 stars. The Mazl5 dataset was updated with 
the KOI data from KOI16. This included the disposition of 
each KOI: ’’false positive”, ’’candidate”, ’’not dispositioned”, 
and ’’confirmed.” KOIs with ’’false positive” disposition were 
removed, as were KOIs with missing rotation periods and 
planet size. These cuts resulted in a list of 2,332 KOIs in 
the Mazl5 dataset. 

From the equatorial coordinates listed in our revised 
McQ14 and Mazl5 datasets, we calculated the Galactic lon¬ 
gitude and latitude coordinates of each target. The Kepler 
held includes target stars spread over a large range in dis¬ 
tance from the Sun, height above the Galactic mid-plane, 
and interstellar extinction; however, they sample only a 
modest range in Galactocentric distance. In order to control 
for differences in Galactic location in our analysis below, we 
calculated the distance of each target using the following 
procedure. 

First, we employed the 2012 version Dartmouth stellar 
evolution isochronesj for a constant age of 4 Gyrs in order 
to derive a simple relationship (in the form of a polyno¬ 
mial) between the absolute K magnitude (Mk) and T e ff and 
[Fe/H]. Our choice of the K-band magnitude for calculat¬ 
ing distance is motivated by the fact that it is the longest 
wavelength magnitude available for each target; there should 
be relatively little interstellar extinction in this band com¬ 
pared to the others. Our choice for the age is close to the 
mean age of simulated stars in the middle of the Kepler field, 
which we generated using the default settings in TRILEGAL 
1. 6^] This relationship is possible because we are restricting 
our analysis to stars classified as main sequence stars. The 
standard deviation of the residual differences between the 
tabulated values of Mk from the isochrones and our sim¬ 
ple interpolation equation is 0.14 magnitudes. Additional 
uncertainties in our estimates of Mk come from uncertain¬ 
ties in the stellar parameters; the typical uncertainty in T e fj 
among the Kepler targets translates into an uncertainty in 
Mk of 0.15 magnitudes. We calculated Mk for each star 
in the McQ14 and Mazl5 datasets, as well as the distance, 
neglecting extinction. 

With a distance estimate in hand for each target, it 
becomes possible to estimate the extinction. The total ex¬ 
tinction along a particular sightline for a specific Galactic 
longitude and latitude was obtained from the NAS A/IPAC 
online calculator j^The Schlafly & Finkbeiner (20111 redden¬ 
ing relation and standard R value (= Ay/E(B — V) = 3.1) 
were used. The ratio of visual to K-band extinction was set 
to 8.8 ( Indebetouw et ah| |2005). Finally, the amount of ex¬ 
tinction in each direction was assumed to reach its maximum 
value 3 kpc from the Sun; larger distances were not required, 
since the target stars are significantly distant from the mid¬ 
plane. The extinction for each target was scaled according to 
its distance relative to 3 kpc. The largest values of K band 
extinction were calculated to be about 0.22 magnitudes, but 


most values were less than 0.1 magnitudes. Nevertheless, we 
corrected the original distance estimates for extinction as 
our final step. 


3 CAMPARING PLANET HOSTS AND 
NON-HOSTS 

3.1 Rotation Period 


In Gonzalez (20111 we compared vsini values of 99 SWPs 


and 594 non-SWPs using primarily the [Valenti and Fischer| 
12005) dataset. The SWPs among nearby stars have been 


detected with the Doppler method. This is in contrast to 
the Kepler field SWPs (classihed as KOIs), which have been 
discovered with the photometric transit method. The nearby 
SWPs have accurate parallax measurements available, while 
very few of the Kepler field targets are close enough for 
parallax distance determination. 

We will apply a modified version of our method of anal¬ 


ysis described in Gonzalez (2008) Gonzalez et al. (20101 and 


Gonzalez (2011) to the present data. In brief, in the prior 


studies we defined an index, called Ai, which is a measure 
of the distance between two stars in T e fi-log g-[Fe/H]-My 
parameter space. This index was calculated for each SWP 
and non-SWP and formed the basis for comparing stars from 
these two groups. 

We now define the following new index, A/+, which is 
more suitable for our present purposes: 

A£ c = 30 | log T e c ff - log T p ff | + | [Fe/Hf - [Fe/H] p \ 

+0.5 | log 5 c - log fl p | + \d c - d p | + \z c - z p | 

where d is the distance from the Sun in kpc, z is the ver¬ 
tical distance from the mid-plane in kpc, p is the KOI index, 
and c is the comparison star (non-KOI) index. If A^ c = 0 
for a given pair of stars, then they are indistinguishable in 
these parameters. 

Our decision to include the parameters d and z is moti¬ 
vated by the need to control for known Galactic-scale trends. 
These include the dependence of the mix of Galactic popula¬ 
tions (thin disk, thick disk, halo) on z, and the concomitant 
trend of metallicity with z. In addition, thick disk stars have 
a different mix of a elements relative to Fe compared to thin 


disk stars; Gonzalez (2009) argued that differences in a ele¬ 


ment abundances should not be neglected when comparing 


SWPs and non-SWPs. For this reason, Gonzalez (2011) em¬ 


ployed [M/H] in their A index definition. According to the 
TRILEGAL simulation of the Kepler field, only about 3 % of 
the main sequence stars should be thick disk stars. There¬ 
fore, their contribution should be almost negligible in our 
analysis. In addition, the Galactic radial metallicity gradi¬ 
ent creates a metallicity trend with d, but it is a small factor 
for our sample. And, as noted above, extinction depends on 
d as well. Although [Fe/H] is explicitly included as a param¬ 
eter in the definition of A^ jC , its uncertainty is large for the 
typical target. However, two stars are more likely to have 
similar values of [Fe/H] if they have similar values of d and 


2 http://stellar.dartmouth.edu/models/grid.html 
2 http://stev.oapd.inaf.it/cgi-bin/trilegal 
4 http://irsa.ipac.caltech.edu/applications/DUST/ 


We employ the Ap iC index in a way similar to that used 


to produce Figure 7 of Gonzalez (2015), which is a compari¬ 


son of lithium abundances in nearby SWPs and non-SWPs. 
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In the present application, the index is used to find which 
non-KOIs in the McQ14 dataset are most similar to each 
KOI in the Mazl5 dataset. The Ap iC index was calculated 
for every pairing of a KOI and a non-KOI and then ranked. 

Next, we classify each KOI star according to the size 
of its planet. If it is between 5 and 15 Earth radii (Re), we 
classify it as a “giant.” If it is between 2 and 5 Re, we call 
it a “Neptune.” We classify smaller planets as “Earths.” 

During the review cycle of the present work, Kepler 
Data Release 24 (DR24) was released. It is the result of 
the first uniform processing of the entire Kepler dataset and 
the first full automation of the dispositioning process. The 
disposition ’’not dispositioned” was not retained in DR24. 
As shown by (Fressin et al. ]2013(|, the rate of false positives 
for giant planets in the Kepler data they analyzed was about 
17%, and it was near 10% for the other planet classes. Given 
the importance of eliminating as many false positives as pos¬ 
sible and of this latest data release, we added one final step 
in our vetting/culling process. We cross-checked our list of 
KOIs in the Mazl5 dataset with the DR24 dataset in order 
to update the dispositions. After removing the ’’false posi¬ 
tive” disposition KOIs for the giant planet class, we are left 
with 108 KOIs with ’’candidate” dispositions and 43 KOIs 
with ’’confirmed” dispositons. We show the results of the 
comparison for the giants in Figure 1. Figure la shows the 
difference in the stellar rotation periods between each KOI 
star and the most similar comparison star. 

Given the relatively large uncertainties in the stellar pa¬ 
rameters for the Kepler targets, we added an additional step 
in the comparison to improve the statistics. Hence, Figure 
lb is the same as Figure la except that the period of each 
KOI star is now compared to the average of the periods of 
the three most similar non-KOI comparison stars. Upon vi¬ 
sual inspection, there appears to be relatively little change 
in Figure lb relative to Figure la, except that the largest 
negative differences are reduced in magnitude in Figure lb. 
It is very clear from both plots that the KOI stars tend to 
have longer rotation periods than the comparison stars. The 
mean differences in Figures la and lb are 17.6 ±2.1 days 
and 17.7 ±2.1 days, respectively]^] If we compare only the 
43 confirmed KOIs, the difference is only slightly reduced 
(19.8 ± 5.2 and 18.9 ± 4.6 days, respectively). 

We show the results for the KOIs we classify as Nep- 
tunes in Figures 2, following the same kind of analyses as 
those of Figures 1. Figure 2 includes 725 KOI stars with 
Neptunes. The mean period differences in Figure 2a and b 
are both 18.0 ± 0.9 days. The mean period differences for 
the 145 confirmed stars with Neptune planets are 18.8 ±2.1 
and 18.4 ± 1.9 days, respectively. 

Finally, the KOIs classified as Earths are shown in Fig¬ 
ure 3. Figure 3 includes 878 stars. The mean period differ¬ 
ences are 20.5 ± 0.9 days and 20.6 ± 0.9 days, respectively. 
The mean period differences for the 157 confirmed stars with 
Earths are 18.8 ± 2.2 and 18.1 ±2.1 days, respectively. 


5 Here and throughout this paper, unless otherwise stated, a 
quoted uncertainty is the standard error of the mean. 



Figure 1. a Difference in the stellar rotation periods between 
each KOI star with a giant planet and the most similar non-KOI 
star, b Same as a but comparing the rotation period of each KOI 
star with a giant planet to the average rotation period of the three 
most similar non-KOI stars. See text for details. 


3.2 vsini 


Before we can compare the results of our analysis of the 
stellar rotation periods of the Kepler field stars to the anal¬ 
ysis of vsini values of nearby stars, we must account for 
the various differences in the way the two samples were cre¬ 
ated. To summarize, homogeneous nearby SWP and non- 
SWP samples, such as those described in |Valenti and Fis^l 
cher (2005 I, are based on planet detections with the Doppler 


method. Doppler surveys are most sensitive to massive plan¬ 
ets with short orbital periods. Essentially all the planets of 


the nearby SWPs sample included in the Gonzalez (2011) 
study would fall under the “giant” planet classification of the 
present work. The non-SWPs comparison sample adopted 


Gonzalez (2011) should contain very few giant planets. 


Therefore, the SWP and non-SWP samples should be rela¬ 
tively pure when the incidence of giant planets is being con¬ 
sidered. When comparing rotation velocities of SWPs and 
non-SWPs among nearby stars in a homogeneous way, the 
best we can do is compare vsini values. 

The Kepler field samples are very different from the 
nearby stars samples. Only a small fraction of planets or¬ 
biting stars in the Kepler field are detectable with the pho¬ 
tometric transit method. The highest fraction of planet de¬ 
tectability occurs for the “hot Jupiters”; about 10 % of stars 
with such planets have detectable transits. This means that 
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Figure 2. Same as Figure 1 but for Neptunes. 


Figure 3. Same as Figure 1 but for Earths. 


any non-KOI comparison sample drawn from the Kepler 
field targets will be contaminated with many SWPs. This 
is simply unavoidable with this type of survey. In the case 
of the giant planets, the situation is not as hopeless as it 
seems, since the statistically inferred true incidence of stars 
hosing giant planets (as we’ve defined them) with orbital 
periods < 50 days in the Kepler field is about 3 % (Howard 
et al.|[2012 l. Thus, something less than 3 % of the stars in 
our non-KOI comparison sample will actually host a giant 
planet. In our McQ14 comparison dataset (with over 15,000 
stars), then, there are nearly 500 stars with undetected giant 
planets. Given this, it is possible that some of the compar¬ 
isons in Figures 1 and 2 are between pairs of SWPs. The 
situation is worse for the Neptunes and Earths, as they are 


much more common; Howard et al. (2012) estimate that 


Neptunes make up 13 %, and Petigura et al. (20131 find a 


similar fraction for Earths. However, our focus here is on 
giant planets. 

In order to compare our results from Figure 1 to the 
results from Gonzalez (2011), we need to convert each stel¬ 


lar rotation period to vsini. Our procedure is as follows. 
First, we adopt a probability density function (PDF) de¬ 
scribing the distribution of the observed inclination angles 
of a population of stellar rotation axes; they are distributed 


like sini (Ho & Turner 2011). We follow Lopez & Jenkins 


(2012) and dehne a PDF for sini using their equation 6 in 


order to produce random samples of sini. We adopt a mini¬ 
mum inclination angle of 5 degrees; SWPs with inclination 
angles smaller than this limit would be difficult to detect 


with the Doppler method. We produced sets of random sini 
values for the McQ14 and Mazl5 datasets and then used 
them to calculate vsini for each star from its rotation pe¬ 
riod and radius. Finally, we added a fictitious error to each 
vsini value drawn fro a Gaussian PDF with a standard devi¬ 
ation of 0.5 km s" , which is the typical uncertainty of vsini 
determinations in the nearby stars samples. Finally, we set 
any resulting vsini values that were less than 0.3 km s" 1 to 
0.3 km s -1 in order to be consistent with the approach of 


Gonzalez (2011). We eliminated vsini values greater than 20 
km s _i , which is necessary because planets are difficult to 
detect with the Doppler method around fast rotators. 

We show the differences in simulated vsini between the 
KOI stars with giant planets and non-KOIs in Figure 4, 
which includes 138 stars. The mean differences in vsini in 
these two datasets are —1.3 ± 0.3 km s -1 . This is just over 
twice the difference quoted in Gonzalez ( 201 1| ) , —0.46±0.11 
km s" 1 , and about twice the value quoted in Gonzalez et 
al. (20101, —0.66 ± 0.13 km s -1 . The agreement is not bad, 


considering the very different natures of the two samples 
and the assumptions that went into producing the Kepler 
simulated vsini dataset. The difference could also be due 
to the fact that the Kepler KOI stars should be older, on 
average, than the nearby stars SWPs, given the Kepler stars’ 
greater distance from the Galactic mid-plane. 
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Figure 4. Data from Figure la transformed to vsini. b See text 
for details. 


Figure 5. Rotation periods of the non-KOI stars (dots) and the 
KOI stars with giant planets from Figure la (plus signs) are plot¬ 
ted against effective temperature. 


4 DISCUSSION 

The results of our analysis of Kepler KOI stars confirms 
previous studies of nearby SWPs indicating that they ro¬ 
tate more slowly than stars without known planets. Our 
results are also consistent with the more narrow finding of 


McQuillan et al. (2013) that close-in planets are rare among 


rapidly rotating host stars. To further explore the nature of 
this difference among the Kepler stars, we show in Figure 5 
the rotation periods of the KOIs with giant planets and the 
non-KOI comparison stars. The first notable difference is ev¬ 
ident among the slowest rotating stars. Stars in the McQ14 
comparison stars dataset have a maximum rotation period 
of 70 days, while the KOI stars with giant planets sample 
have periods up to 142 days; there are 16 KOI stars with 
giant planets with periods longer than 70 days. Just above 
the ridge of greatest concentration of stars, the KOI stars 
have a high relative incidence. In addition, among the fastest 
rotators, KOI stars have a very low relative incidence. 

The McQ14 and Mazl5 datasets both make use of the 
same analysis methods on the same Kepler data (Q3-Q14). 
Thus, it is unlikely that the differences between the KOI 
and non-KOI stars evident in Figure 5 are due to differences 
in the data analysis. Nevertheless, as a test of the reality 
of the difference between the mean KOI and non-KOI stars 
rotation periods, we have calculated a new mean period dif¬ 
ference for the data in Figure la by excluding the 16 KOI 
stars with giant planets with periods greater than 70 days; 


we obtain 11.4 ± 1.6 days. This is still a highly significant 
difference. 

The present results are particularly relevant to research 
on the lithium abundances in SWPs, which has received 


more attention (Figueira et al. 2014 Gonzalez 2015 Del 


gado Mena et al. 20151. citetis04 suggested that since the 
protoplanetary disc around a young star contains a large 
fraction of the systems angular momentum, it can cause ro¬ 
tational breaking of the star. The breaking, in turn, results 
in deeper mixing in the star and thus more efficient destruc¬ 


tion of Li. Bouvier (2008) placed this idea on a firmer phys¬ 


ical foundation. The implications of this idea are important 
for both planet formation and stellar evolution (especially 
gyrochronology) theories. If correct, Bouvier (20081 would 
place a planet-hosting star’s spin-down very early in its his¬ 
tory, while it still possessed a protoplanetary disk. 

There is strong observational support for the slowing of 
stellar rotation with increasing age for main sequence stars 
younger than about 1 Gyr ( |Soderblom|20ld) |Meibom et al.| 
20151. Among main sequence stars in clusters with ages less 


than 300 Myr, the measured dispersion in rotation periods 
is very large, about 0.1 to 10 days. By the age of the Hyades 
and Praesepe (~ 600 Myr), the observed dispersion in rota¬ 
tion period is much smaller (Del orme et al.|2 011). 

Observational evidence for stellar slowdown among 
older stars is sparse, given that rotational modulation due to 
the presence of active regions becomes weaker with increas- 
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ing age. This, again, is where Kepler data are particularly 
valuable. Meibom et al. (2015) have measured the rotation 
periods of 30 faint (all fainter than V = 15 magnitude) main 
sequence stars in the 2.5 Gyr old cluster, NGC 6819. Theirs 
is an important study of stellar rotation, as it extends di¬ 
rect measurement of stellar rotation periods to ages much 
greater than 600 Myr. Their data show a steady trend of 
increasing rotation period as one goes down the main se¬ 
quence with relatively small scatter. However, there is still 
some scatter in the periods around a mean trend curve that 
is greater than the quoted uncertainties in the rotation pe¬ 
riods (see their Figures 2 and 3). In particular, the Kepler 
stars 5112268 and 5025271 have nearly identical colors and 
magnitudes and yet differ in rotation period by nearly 3 
days (with period uncertainties of 0.10 and 0.25 days, re¬ 
spectively) . 

If our main result is correct, then about 30% of the 
main sequence stars in NGC 6819 should be accompanied 
by planets and thus have longer rotation periods than stars 
without planets. This should result in a large scatter of the 
periods measured by Meibom et al. (20151. At first glance, 
the relatively small scatter found by them would seem to 
be inconsistent with our findings. We can offer two possible 
explanations for the apparent inconsistency. First, there is 
an observational selection bias against detecting stars with 
longer rotation periods. Slowly rotating main sequence stars 
are chromospherically less active than faster rotating stars, 
which makes it more difficult to measure rotation from star 
spots rotating into and out of view. A possible way to test 
this idea is to calculate the variability amplitude of main 
sequence stars with rotation periods 10 to 50 days longer 
and estimate what fraction would have been detectable in 
their analysis. 

Second, the environment in a rich cluster like NGC 6819 
might not be representative of the kind of birth environment 
of the typical star in the Kepler field. Only a handful of 
planets have been found in clusters. We need better statistics 
to test this idea. 


5 CONCLUSIONS 


We have verified that the rotation periods of KOI stars in the 
Kepler field differ significantly from those of otherwise sim¬ 
ilar non-KOI stars. Stellar rotation period measured with 
precision photometry is preferable to spectroscopic vsini 
measurements when doing this kind of comparison, given 
the small relative errors that can be achieved. Neverthe¬ 
less, we have also shown that the rotation period measure¬ 
ments of KOI stars with giant planets are consistent with 
nearby SWP vsini measurements. Both samples show that 
stars with giant planets rotate more slowly than stars lack¬ 
ing such planets. 


Bouvier (20081 suggested that slow rotation among 


SWPs caused by early star-disk interactions caused accel¬ 
erated Li destruction in them. Recently, additional evidence 
for lower Li abundances among SWPs has been presented 


by Figueira et al. (2014 ) and Gonzalez (20151. These results, 


as well as the new findings presented in this work, further 
strengthen the links among Li, rotation and the presence of 
planets. However, additonal study is required to determine 
the timing of the spin-down of the host stars. Bouvier’s the¬ 


ory requires that the spin-down would occur mostly very 
early in a star’s history. What is surprising is that compa¬ 
rably slow rotation is also found among stars in the Kepler 
field hosting Neptunes and Earths. Also myterious within 
the context of these findings is the apparently small scatter 
in measured rotation periods of main sequence stars in old 
clusters. 

The next steps in this line of research could involve 
correction for the presence of stars with planets (especially 
Neptunes and Earths) in the comparison samples. This kind 
of contamination is unavoidable in photometric transit sur¬ 
veys, but it should be possible to account for it in a statisti¬ 
cal way. In addition, more rotation period measurements in 
star clusters are needed in order to determine if the scatter 
in rotation period among old stars is really as small as it 
appears to be from current data. 
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